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In the present work, it is developed a formalism to deal 
with the macroscopic study of the astrophysical systems, 
which is based on the consideration of the exponential self- 
similarity scaling laws that these systems exhibit during the 
realization of the thermodynamic limit. Due to their scaling 
laws, these systems are pseudoextensive, since although they 
are nonextensive in the usual sense, they can be studied by 
the Boltzmann-Gibbs Statistics if an appropriate representa- 
tion of the integrals of motion of the macroscopic description 
is chosen. As example of application, it is analyzed the sys- 
tem of classical identical particles interacting via Newtonian 
interaction. A renormalization procedure is used in order to 
perform a well-defined macroscopic description of this sys- 
tem in quasi-stationary states, since it can not be in a real 
thermodynamic equilibrium. Our analysis showed that the 
astrophysical systems exhibit self-similarity under the follow- 
ing thermodynamic limit: E ^ oo, L ^ 0, N ^ oo, keeping 
E/N3 =const, LN3 =const, where L is the characteristic 
linear dimension of the system. It is discussed the effect of 
these scaling laws in the dynamical properties of the system. 
In a general way, our solution exhibits the same features of 
the Antonov problem: the existence of the gravitational col- 
lapse at low energies as well as a region with a negative heat 
capacity. 

PACS numbers: 05.20.-y; 05.70.-a 



I. INTRODUCTION 

Traditional thermodynamics is not able to describe the 
astrophysical systems. They do not fulfill the additivity 
and homogeneity conditions, indispensable requirements 
for the good performance of this formulation. 

In the last years it has been devoted so much effort to 
the extension of the Thermodynamics to the study of the 



nonextensive systems. In this context, the astrophysical 
systems have received a special attention. 

Among the remarkable new results obtained in the 
present frame, it can be mentioned the application of 
the so popular Tsallis' nonextensive statistics to the 
analysis of astrophysical systems. In this approach it has 
been justified the application of the polytropic models 
1^,^, which have been extensively used in the descrip- 
tions of such a systems H,p| : 



p = Cp^ ,j>l, 



(1) 



where p is the pressure, p, the particles density, and 7, 
the polytropic index. 

However, in our opinion, in those works it has been put 
off the microscopic justification of the polytropic models 
to the applicability of the Tsallis' theory. The nonex- 
tensive statistics is not a completely satisfactory formu- 
lation, due to the all theory dependence on the entropic 
index, q, which is a parameter representing a measure of 
the degree of nonextensivity of the system. In spite of 
the attractiveness of this formulation, the theory is not 
able to determine univocally the value of the entropic 
index, at least, in the context of the equilibrium Ther- 
modynamics, so that, it must be appealed to the exper- 
iment or computational simulations in order to precise 
it. Some evidences aim that the entropy index could be 
obtained throughout the sensitivity of the system to the 
initial conditions and the relaxation properties towards 
equilibrium 

Alternative approaches have been proposed using the 
microcanonical ensemble. Although it canjiot be assured 
its application to any Hamiltonian systemEJ, the Thermo- 
statistics could be justified starting from this ensemble 
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^Apparentely the microcanonical ensemble is only well- 
defined statistical ensemble, whose justification ordinarily is 
based on the chaotic properties of the trayectories for a generic 
non-integrable system when it is overcome a few tens of de- 
grees of freedom. However, there are some computational ev- 
idences in dynamical studies of some nonextensive systems 
in which the mixing time is extremely long, and diverges 
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without invoking anything outside the Mechanics. This 
description is apphcable to many situations in which the 
canonical description fails, allowing us to determine the 
necessary conditions for the applicability of any general- 
ized canonical ensemble. 

Starting from this viewpoint, in the ref. it was ad- 
dressed a generalization of the extensive postulates of the 
traditional Thermodynamics in order to extend its ap- 
plication to the nonextensive systems. According to our 
proposition, this objective could be carried out taking 
into consideration the self- similarity scaling postulates: 
the equivalence of the microcanonical ensemble with a 
generalized canonical one during the realization of the 
thermodynamic limit throughout of the self-similarity 
scaling properties of the system fundamental physical 
observables: the behavior of integrals of motion, the ex- 
ternal parameters and the accessible volume of the micro- 
canonical description with the increasing of the system 
degrees of freedom. 

So far, these postulates have been applied to the anal- 
ysis of the necessary conditions for the validity of two 
statistical formulations: in the ref. p^ , to the micro- 
canonical thermostatistics of D. H. E. Gross [ p^UTs] ], as 
well as in the ref. to the Tsallis' nonextensive statis- 
tics |. 

The first is a theory based on the consideration of 
the microcanonical ensemble with the assumption of the 
Boltzmann's definition of entropy: 

SB^\nW, (2) 

his famous gravestone epithaph in Vienna. Since the 
Boltzmann's entropy does not demand the realization of 
the thermodynamic limit this formulation is applicable 
to some small and mesoscopic systems. The thermody- 
namic formalism of this theory has been defined in order 
to be equivalent to the traditional one when it is ap- 
plied to extensive systems. That is the reason why this 
theory is appropriate to the macroscopic description of 
those systems becoming extensive when the thermody- 
namic limit is invoked although they are not found in 
the thermodynamic limit. 



with the increasing of system degrees of freedom. In some 
of these cases is also observed the existence of anomalous 
quasi-stationary states which can not be described using the 
microcanonical ensemble. Only after a long trasient time, the 
system abandons these singular states and starts slowly ap- 
proaching to the equilibrium dictated by the microcanonical 
ensemble. In this case it is shown that the infinite time limit 
{t — > oo) necessary to the stablishment of the equilibrium de- 
scribed with microcanonical ensemble does not commute with 
the thermodynamic limit (A'^ — > oo), since the second is per- 
formed before the first, the quasi-stationary state becomes the 
real equilibrium state of the system. The interested reader 
may check the following reference for details: Phys. Rev. E 
64, 056134 (2001). 



In the ref. jis) it was shown that the Gross's theory is 
also applicable to all those Hamiltonian systems exhibit- 
ing exponential self- similarity scaling laws in the thermo- 
dynamic limit. We called these system as pseudoexten- 
sive when N —f oo, since in spite of they are nonexten- 
sive in the usual sense, they can be tried by means of 
the Boltzmann-Gibbs' statistics if an appropriate repre- 
sentation of the integrals of motion is chosen. Many of 
the systems found in the real world belong to the class 
of the pseudoextensive systems, since it is enough an ad- 
ditive kinetic part in the Hamiltonian of the system for 
exhibiting this kind of scaling laws in the thermodynamic 
limit. 

In the ref. was shown that the Tsallis' statistics 
is appropriate to the macroscopic description of those 
Hamiltonian nonextensive systems exhibiting a potential 
self- similarity scaling laws in the thermodynamic limit. 
In our approach many details of this formalism natu- 
rally appear, starting only from the Mechanics under 
the consideration of this kind of scaling laws, i.e., the 
g^expectation values, the ^-generalization of the Legen- 
dre's Transformation (see for example in refs. p^ , p^ ). 

The above reasons allow us to consider that the astro- 
physical systems belong to the class of the pseudoexten- 
sive systems, and therefore, it is also justified the appli- 
cation of the Boltzmann-Gibbs' statistics to the study of 
these systems, at least, when the phase transitions are 
not present. 

The Tsallis' statistics is expected to describe systems 
exhibiting potential distributions, that is, systems with 
fractal characteristics. In the ref. it is suggested that 
the Tsallis' potential distribution fitted very well the dif- 
ferential energy distribution of the dark matter in halos 
obtained by a numerical simulation. Similar analysis have 
been carried out in ref. by Fa & Pedron for the ellipti- 
cal galaxies. However, the fractals properties can be also 
found in the context of Boltzmann-Gibbs' statistics, for 
example, in the Michie-King models for globular clusters 
1^ (see ref. ||2|] for review), or the models proposed by 
Stiavelli & Berlin ||2|LHjorth & Madsen (24) for elliptical 
galaxies, and others [p5[ . These models lead to composite 
configurations with an isothermal core and a polytropic 
envelope and can not be justified by the Tsallis' gener- 
alized Thermodynamics. Many authors state that these 
models are more appropriate for the description of those 
astrophysical objects (see for example in ref. p6[|). 

In the present work it is pretended to reconsider the 
thermo-statistical description of the astrophysical sys- 
tem, but this time, taking into consideration the self- 
similarity scaling postulates in order to find the 
necessary conditions for the validity of the generalized 
canonical ensemble (Boltzmann-Gibbs) in the thermody- 
namic limit. 
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II. MICROCANONICAL DESCRIPTION. 
A. Microcanonical Mean Field Approximation. 

Let 5 be a Hamiltonian system composed by a huge 
number N of identical particles, which interact among 
them by means of short-range and long-range forces si- 
multaneously. Let us also consider that the character- 
istic linear dimension of the system is comparable with 
the effective radio of the long-range interactions, but it 
is extremely large in comparison with the effective radio 
of the short-range interactions. 

The above conditions allow us to speak about two small 
scales in the system: 

• Microscopic scale: linear dimensions comparable 
with the effective radio of the short-range interac- 
tions. 

• Local scale: very large linear dimensions in regard 
to the microscopic scale, but extremely short in 
comparison with the effective radio of the long- 
range interactions. 

Let us consider a partition of accesible physical space 
in not overlapped cells, {cfc}, whose characteristic linear 
dimensions correspond to a local scale. It is easy to show 
that the iV-body phase space integration can be descom- 
posed in the following way: 



1 

m 



dX 



N 



N 

k rifc=0 



5^^'^ N- 




(3) 



where it have been taken into account all the possible 
configurations for distributing N identical particles in the 
cells. In this expression, Xm represents the A^-body phase 
space, while A'i^'' represents n^-body phase space whose 
accesible physical space have been limited to the k-th 
cell, Cfc. In addition, it has been introduced the following 

integral operator O 



O 



i / dX„, if n^O, 



1, if n = 0, 



(4) 



as well as the function 5'^'^' [n] = 6on, which assures the 
particles number conservation. A derivation of the Eq. (|^) 
appears in the appendix A. 

The enormous scale separation among the effects of 
the short and long-range interactions supports the valid- 
ity of the spacial adiabatic approximation. Let us con- 
sider those configurations in which each cell contains 
a very large number of particles n^. and let us denote 
this subsystem of particles by Sk- The physical quan- 
tities characterizing the long-range interactions, that is, 
the long-range interacting fields (r), almost do not vary 



at the spacial region occupied by the cell due to its lin- 
ear dimension. It means that the long-range interactions 
almost do not distinguish the internal structure of the 
subsystem enclosed by the cell, so that, these interac- 
tions are only effective for the subsystem as a whole. 

Therefore, the contribution of long-range interactions 
to the system total energy can be approximated by func- 
tional terms of the mean values of the long-range inter- 
acting fields at the region of each cell, (f>k, as well as 
of certain collective quantities qk characterizing the sub- 
systems {iSfc}, like the mean values of particles density, 
magnetization density, etc. Similarly, the mean values 
of the long-range interacting fields can be determined 
from the mean values of some collective quantities of 
the subsystems {Sk}. On the other hand, due also to 
the linear dimensions of the cells, it can be neglated the 
short-range interactions among the particles beloging to 
different cells. 

Taking into account all the above exposed, the system 
Hamiltonian can be approximately expressed as follows: 



H ^int (^«.; + Vioc (h, qk) , 



(5) 



where (^fe; (t>k) is the internal energy of the n^-body 
subsystem Sk enclosed in the cell Ck , which only involves 
the kinetic energy of the particles as well as the contribu- 
tion of the short-range interactions among them, being 
Xn^ their microscopic degrees of freedom. Here, it is also 
included a parametric dependence of the internal energy 
of the fields in order to take into account a possible 
influence of these fields on the internal configurations of 
the subsystem Sk- On the other hand, Vioc is a local 
term of energy containing the contribution of the sub- 
system Sk as a whole, which only involves its effective 
interactions with the long-range interacting fields (jfk- As 
already mentioned, the fields 0fc are determined through 
a determined functionals of certain set local quantities, 
q = {qk}, which characterize the local subsystems: 



(l^k = J^k {q) ■ 



(6) 



As it could be seen, the effective long-range interact- 
ing fields 0fe at the k-th cell can be considered as external 
parameters for the subsystem Sk. Therefore, each sub- 
system Sk can be considered as locally extensive. Due to 
the presence of the long-range interactions, the system 
5 = y 5fe is nonextensive: the quantities characterizing 

the subsystems as a whole vary during the continuous 
passage among the cells. 

There are systems in which the energy is not the only 
one integral of motion determining their macroscopic de- 
scription (for example: the astrophysical systems), so 
that, it could be considered other integrals of motion, 
like total angular momentum ||2l[]. In such as cases, it is 
well-known that the total angular momentum admits a 
descomposition similar to the expression of the energy in 
the Eq.(|): 



3 



M 



+ Ukrk X Vfc, 



(7) 



where and are the position and velocity of the mass 
center of the subsystem Sk, fJ^k is its total mass, while 

•^irat^ and Xn^. are the internal angular momen- 

tum and the microscopic degrees of freedom of the sub- 
system in its own mass center frame. The consideration 
of the collective motion of the subsystem Sk leads to add 
a kinetic term ^l^k^l to the sum in the Eq.(|^) as well as 
to substitute by . 

Thus, the integrals of motion involved in the macro- 
scopic description of the system can be represented in the 
following form after the considerations assumed above: 

/ -X! ^internal i^n^ 'Ak) + Jcollective ((/>fc , q) , (8) 
k 

where ^l^^t^^^ai i^n^i 4'k) represents the internal contri- 
bution from the microscopic degrees of freedom of the 
subsystem Sk, while Jcoiiective{4'k,qk) is the contribu- 
tion of its collective degrees of freedom at local level. 

According with what was exposed above, the accessible 
volume of the system in the microcanonical ensemble is 
approximately given by: 

W (I, N) d^ki E a [Xi'^] S - Tk (q)] X 

k \nfc=0 



k 



nk 



xS 



I- E ^internal (^Kfc ' <^fc) + Jcollective (0fc, qk,^k) 



(9) 

where it has been explicitly introduced the dependence of 
the collective term of the mass center velocity, in order to 
take into account the collective motion of the subsystems. 
Using the identity: 



dikdvkS 



Ik 



ai (-'^"fc ; '^fc) 









)1 


s 


Vfc - — Pfe 









1, 



(10) 



where Pk is the total linear momentum of the sub- 
system Sk, the Eq.(^ can be rewritten introducing 
the Boltzmann's entropy for each local subsystems Sk, 

SB{ik,nk;<l>k)- 

exp [5*5 (ife, rife; = 



^k-l'^Lal{Xn,;ct^k)]5[Vk], (11) 



where ik is the internal contribution of the subsystem to 
the total integrals of motion of the system S. Thus, the 
accessible volume of the system can be computed from 

N 

W{I,N) E / dikdc^kd^rk5[(j>k-H{q)]y- 

k nfc=0 •' 



X exp 



i'k,nk; <Pk) 



L k 



(^N- ^ Uk^ 



collective \H^k i Qk ? j 



(12) 



where the physical quantities qk are determined from the 
ik by means of certain functional dependencies: 



qk = fq {ik,nk,^rk) 
Developing the continuum limit: 



N ^ CO, and ^ [cj) / {^Z ^ ^'^^^ ~^ ^' 



(13) 



(14) 



where n (ck) is the physical volume of the k-th cell, the 
subsystem Sk is locally perceived as a fluid. Thus, it is 
obtained finally the microcanonical mean field approxi- 
mation (MFA).- 

W {I, N) ~ Wmfa {I, N) = 

C jvg (r) Vp (r) (r) Vw (r) 5 {0 (r) - [r; g, p, v] } x 

J d'r s[g{r),p{r);<f>{r)] S (^N - Jd\p{r)^x 

I- y^^r e(r)+,9[r;(/.(r),e(r),p(r),v(r)] , 

(15) 



X exp 



x^ 



where p(r), g{r), i9 [r,0 (r) , p (r) , p (r) , v (r)], and 
s [g (r) , p (r) ; (/) (r)] are respectively: the particles den- 
sity, the densities of the internal and collective contri- 
butions to the total integrals of motions of the system, 
and the entropy density of the fluid at the neighborhood 
of the point r. C is an unimportant constant which ap- 
pears as consequence of the continuous limit and can be 
ignored. 



B. The case of the astrophysical systems. 

Let us apply the microcanonical mean field approxima- 
tion to the analysis of an astrophysical system composed 
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by a static fluid of identical particles, whose macroscopic 
state is only determined by its total mass (number of par- 
ticles) and the energy. Let e (r), p (r) and s [e, p; 0] be re- 
spectively the internal energy, particle, and entropy den- 
sities of the system neighborhood of the point r, where (p 
is the sclf-gravitating Newtonian potential. In this case 
the microcanonical mean field approximation is written 
as follows: 

Wmfa {E, N)^ jve (r) Vp (r) V4> (r) 5 {</> (r) - C [r; p]} 



X exp 



j (fir s [e (r) , p (r) ; (r)] S N- J d^r p{r 



xS 



E- j dhH[e (r),p(r),0(r)] 



(16) 



where H [e (r) , p (r) , (/) (r)] is the Hamiltonian density of 
the system, which is given by: 



H [e, p, I 



1 2 

(V0) +mp(j) + e. 



8ttG 



(17) 



In the above definition: G is the Newton's constant and 
m is the particles mass. Finally, £ [r; p] is the functional: 



C[r; p] — m / d r' g (r, r') p (r') 



(18) 



which determines the spacial configuration of the New- 
tonian potential (j) (r) at a given spacial configuration of 
the particles density, /9(r), where (7(r, r') is the Green's 
function of the Poisson's problem: 



AG (r, r') = AttGS (r - r') 
In the tridimensional case: 



(19) 



(20) 



Using the Fourier's representation of the Dirac's delta 
function: 

/dk 
— exp [z (a; - x')] , (21) 
27r 

where z = P + ik, with /? G R, the Eq.(p^) can be rewrit- 
ten as follows: 



Wmfa [E, N) 



dk 
(2^ 



■Z[zi,N]e^p[ziE]. (22) 



The functional Z [zi, N], with argument zi = (3 + ik, is 
given by: 



+ 00 

f ds 

Z [zi, TV] = / — exp [z2N] H (zi, za) 
J 27r 



(23) 



with Z2 = -F is , /.t G R, and the function H (zi, Z2) is 
expressed as follows: 



H(Z1,Z2) 



J Ve (r) Vp (r) Vtj) (r) VJ (r) x 



X exp[-i? (£,p,0,w;zi,Z2)] . (24) 
The functional H {e, p,4>,uj; zi, Z2) is expressed by: 

H {e,p,(j),uj;zi,Z2) = J d^r zin[e (r) , p (r) , (r)] -|- 

+Z2P (r) - w (r) {(f) (r) - £ [r, p]} - s [e (r) , p (r) ; (r)] , 

(25) 

where lo {r) = j (r) -I- iJ{r). The auxiliary field J (r) 
allows us the Fourier's representation of the delta func- 
tional of the Newtonian potential, 0, and j (r) is an ar- 
bitrary real function. 

As it was previously pointed out, the astrophysical 
systems can be considered as pseudoextensive |13|: they 
exhibit an exponential self-similarity scaling laws in the 
thermodynamic limit |p^ , |l3| , the limit of many particles. 
As consequence of this behavior, they can be dealt with 
the usual Boltzmann-Gibbs' Statistics, only if an appro- 
priate representation of the integrals of motion determin- 
ing their macroscopic state is chosen and if they do not 
present first-order phase transitions. In order to select 
a correct representation for the integrals of motion, the 
self-similarity scaling laws of the system must be found. 

It is demanded the following symmetry: 



N aN, p - 
PE a(3E, V 



> aV 



(26) 



in order to make equivalent the microcanonical with the 
canonical description when N 00, where V is the 
Planck's potential: 



r{/3,N) = -\nZ[P,N]. 



(27) 



and Sb, the Boltzmann's entropy, Eq.(g). The above re- 
quests characterize the system as a whole. However, the 
description of the system is performed at a local level, and 
therefore, the scaling laws must be determined for the lo- 
cal fields: £ (r), p (r), </> (r), j (r), s [e, p; (j)], the parameter 
f3, as well as the scaling law of the spacial coordinate r. 
Let us consider the following scaling laws: 



■ ^J■, j (r) J (r) , p^ a'^p. So a™so, 



(28) 



where c, tt, ^ ^^'^ ^-re certain real scaling expo- 
nent constants. From the analysis of the Eq.(p5|) it is de- 
manded the following relations in order to satisfy the ho- 
mogeneous scaling of the functional H (e, p, 0, lo; zi, Z2): 
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3c + m = 1, 3c + >f = 1, 3c + TT + X = 1, 

T^~2c^ K,2{j^-c) = x- (29) 
The solution of the above equation system is given by: 
m = >c, 7]^ TT = -r], 

c=^, X=^- (30) 
The scahng law for the total energy is given by: 

E^a'E, (31) 

where: 

>ir+ 5 



and therefore, an appropriate selection of the represen- 
tation of the integrals of motion is: 



/ = {£,N) with 8 = E/N"^. 



(33) 



In this case, the correct Legendre's transformation be- 
tween the thermodynamic potential is given by: 



with: 



Sb{.£,N)~(3o£-V{13o.N), 



l3o = f3N'^ and£ = —V{Po,N). 



(34) 



(35) 



The above scaling laws depend on the parameter 
so that, they are only specified when the microscopic 
model for the fluid is assumed. When the thermody- 
namic limit is performed in the generalized canonical en- 
semble it will found that the system will carry out more 
probably those configurations minimizing the functional 
H {e, p,4>,uj; zi, Z2). Thus, it is arrived to an equilibrium 
mean field theory. These configurations are obtained by 
solving the following equations: 



6H (£,p, j;/3,^) 
Se (r) 



= 0, 



dpir) 



SH {e,p,(j),j](3,p) 
H (r) 

imposing the contrains: 

SH {e,p,(l),j;P,p) SH {£,p,(j>,j;(3,p) 



0, 



(36) 



Sj (r) 



0, 



Sp 



N, (37) 



which are related with the conservation of the particles 
number and the consideration of the Poisson's equation 
for the Newtonian potential 0. The maximization takes 
place when it is guarantied the positive definition of the 
functional matrix: 



A,(r',r) 



S^H 



Sh{v')Sf,{v) 



(38) 



with /i=i,2 = (£)P)i where the subindex s represents the 
solution for the minimization conditions, Eqs.(^6|) and 
(^). From the conditions given in the Eq.(36) are de- 
rived the following relations: 



d 



(39) 



M= (Ss, Ps] (ks) - f3m(j)s - C {r, j,) , (40) 

OPs 

the conditions for thermodynamic equilibrium along the 
volume of the system, as well as the structure equation: 



(32) A(j)s = AttG 



d 

mps - I3~^ ( js + TTT* {^s,Ps\ 4>s) 



(41) 



From the validity of the Poisson's equation: 

A(t)s ^ A-KGmps, (42) 

the functional dependency for the auxiliary field jg (r) is 
derived: 



3s 



d 



s [ss^Ps; 4>s) 



(43) 



It could be also used another alternative choice for the 
Hamiltonian of the system given in the Eq.(p^), : 



Ti. (e, p, ( 



e + —mpcf). 



(44) 



In this case equilibrium conditions lead to the following 
relations: 



d 

P = ■Q^s{es,Ps;(j)s) , 

d 1 

M = {es,Ps;(ps) ~ ^fim(l)s - C{r,js) 
ops 2 



1 



-dm(f)s 



d 



s {es,Ps; (t>s) 



(45) 
(46) 
(47) 



Taking into consideration the Green's solution for the 
Poisson's equation, the Eq.([l^), it is deduced the rela- 
tion: 

d f 9 \ 

= -g^S {Es, PslCps) - Pm(j)s - C ( r, --g^S {Es, Ps] (f>s)j , 

(48) 

which is the same obtained as using the first Hamilto- 
nian, the Eq.([2^). Let us introduce the function C (r) as 
follows: 



C(r) -£ r 



d 



s ies,Ps;4>s) 



(49) 
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Summarizing: the equations that dictate the equihbrium 
of the system are the following: 

d d 
P = ^— s {es,Ps;4's) , M = {es,Ps;(t)) - I3m(j)s + C, 

(50) 



A(/., = 47rGmp„ AC = 47rG— s (e,, p,; 0,) . (51) 

0(ps 

which have to satisfy the minimum request: the non- 
negativity of the matrix functional of the Eq. (|3^) . It 
is also necessary the exigency of the conservation of the 
number of particles: 

j cPrps (r) . (52) 

The Planck's potential is expressed in terms of the local 
description as follows: 

J d^r /3H [Es (r) , (r) , 0, (r)] - s [e, (r) , (r) ; 0, (r)] , 

(53) 

where /3 is related with 0o by the first relation of the 
Eq. (^5|) . Using the Eq.(^J), the above relation is rewrit- 
ten introducing the function p [f3, p\ 0]: 

V (/3o, N)^ j dPv p^jip, (r) 0, (r) -f p p, (r) ; 0, (r)] , 

(54) 

which is the Planck's potential density at the point r. 
Due to the scaling laws, it is convenient to set = 1, 
and let the dependency on N to the scaling parameter a. 
That is to consider N as scaling parameter: 

a = TV, (55) 

and set iV = 1 in all the above relations. In this case 
it will be performed a scaling invariant description of 
the system. The equation system, the Eq.(|l]), must be 
solved under the constrain: 

y d^r p, (r) = 1. (56) 

Finally, the energy scaling invariant is given by: 

e (/3) = y" d^r E, (r) + ^mp, (r) 0, (r) . (57) 

In this way, when 'N dependency is taken into account 
in the scaling laws, the Eq. (|28|) , it is obtained a suitable 
problem that could be solved by means of numerical com- 
putation. The microscopic model for the fluid must be 
specified in order to continue the calculation. As example 



of application, let us to apply this formalism in the next 
section to the analysis of the classical system of identical 
non-interacting particles (at the microscopic level) . 

As it can be seen, the consideration of the self- 
similarity scaling postulates do not altere the descrip- 
tion in the scaling invariant equilibrium mean field theory 
{N = 1). However, if the scaling laws are not correctly 
chosen, this fact must lead to some unphysical conse- 
quences, such as the non-proportionality of the Boltz- 
mann's entropy and the Planck's potential with the par- 
ticle number, similarly to the Gibbs' paradox, and there- 
fore, it there will be a trivial ensemble inequivalence. 

III. IDEAL GAS OF PARTICLES. 
RENORMALIZATION. 

The main difficulty in the statistical description of the 
astrophysical systems is the existence of both, a short- 
range and long-range singularities due to the considera- 
tion of the Newtonian gravitational interaction. The first 
situation is the very-well known gravothermal catastrophe 
of the TV-body self-gravitating system [ pTp^ ] . In such a 
system, there is no upper bound on the entropy and a 
state of arbitrarily large entropy can be constructed from 
a centrally concentrated density profile by shifting more 
of the mass towards the center (core-halo structures) . It 
can be seen in p^ , ^ , ^ for review. This situation can 
be easily avoided, since a new Physics appears at micro- 
scopic scales, i.e., the Quantum Physics, which consti- 
tutes a natural renormalization when the system is con- 
stituted by micro-particles: molecules, atoms and sub- 
atomic particles, or in general way, by the consideration 
of the particles size. 

The second, the long-range singularity, has a differ- 
ent nature. It is very well-known that the gravitation 
is not able to confine the particles: it is always possi- 
ble that some of them have the sufficient energy for es- 
caping out from the system, so that, the system always 
undergoes an evaporation process. Therefore, the astro- 
physical systems will never be in thermodynamic equi- 
librium. However, there are intermediate stages where 
this behavior might be neglected and a quasi-equilibrium 
state might be reached (dynamical issues like ergodic- 
ity, mixing or "approach to equilibrium" [ ^2[j3^j3^ ). In 
principle, these quasi-stationary stages can be described 
dismissing the system evaporation. This last process 
could be considered as a secondary effect, which mod- 
ifies the quasi-stationary equilibrium of the system. In 
the present approach, it will be only analyzed the quasi- 
stationary equilibrium. 

Let K be the kinetic energy of a given particle, and 
(j) (r) its correspondent Newtonian potential energy at the 
point r. That particle will be retained by the system 
gravity if the following condition is hold: 

K + m4><0. (58) 
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When the above condition is not satisfied, the particle 
will be able to scape out from the system if it does not lose 
its excessive energy. It will be only consider in our de- 
scription those quasi- stationary stages in which the above 
condition is hold for each particle of the system. This as- 
sumption is the key of the so called Michie-King models 
(see refs. pl| , p^ , p3| , p4| , p5| ) . This exigency acts as a reg- 
ularization procedure, since it is sufficient to confine the 
system, avoiding in this way the long-range singularity. 
No rigid boundaries are necessary in this case, so that, no 
artificial parameters like the box volume are introduced 



The local entropy density, that is, the function 
So (EjP; (t>), is obtained from the following model: Let us 
consider a classical system of TV non-interacting parti- 
cles which is confined by means of a rigid boundary with 
volume V. Taking into consideration the above renormal- 
ization prescription, the Eq.(^8|), the admissible stages of 
this system are those in which the kinetic energy of each 
particle satisfies the condition: 



2m 



(59) 



where m is the particle mass and U, the energy thresh- 
old. The renormalized partition function of the canonical 
ensemble, Z-ji (A, N; U), is given by: 
Zn{X,N;U)^ 



1 d^^Qd^^P 
^ (27rfi)3^ 



(60) 

where A is the canonical parameter. The calculation 
yields: 



P{X,N;V,U) = 



{xuy- 



and therefore, the energy is: 
E 



— (l-l zd,lnF(z)\ 
2A V 3 ^ ^lz=(A(7)5 



(65) 



(66) 



The caloric curve is shown in the FIG. 2. It can be 
seen the linear behavior of the energy, at low values of 
the 'temperature', T — X~^, which corresponds with the 
usual Maxwell's distribution. With the increasing of the 
parameter T, the divergency between the renormalized 
model with the ideal gas system becomes evident. This 
asymptotic dependency characterizes a uniform distribu- 
tion function for the particles momenta. In the FIG. 3., 
it is shown the behavior of the distribution function of 
p for different values of the parameter z = {XUY ■ This 
graphic shows the transition from a gaussian distribution 
for high values of z, to the uniform distribution at value 
z = 0. Finally, the entropy function is obtained by means 
of the Legendre's transformation: 



S{E, N; V,U)^XE~P (A, N; V, U) . 



(67) 



FIG. 1. Behavior of the function F (z). With the increasing 
of z this function tends lastly to the unity. 



FIG. 2. Caloric curve of the microscopic model. At low 
"temperatures", T = A~^, the model behaves like the ordi- 
nary ideal gas, with a maxwellian velocity distribution func- 
tion, but at high temperatures, the energy cutoff acts becom- 
ing homogeneous its velocity distribution function. 



Z^iX,N;U)='-[^^y {F[iXUY]y. (61) 
The function F (z) in the above expression is defined by: 

z 

F{z) = ^ [ exp (-x^) dx, (62) 



that is: 



F (z) = erf (z) exp (— z^) 



(63) 



In the FIG. 1., it is shown the behavior of this func- 
tion. The asymptotic dependency for low values of its 
argument is given by: 



F{z) 



30F 50F 



-f O (z^) , for z < 0.5. (64) 



The Planck's potential in the thermodynamic limit is 
given by: 



FIG. 3. Velocity distribution function at different values of 
z. All these functions were normalized as the unity at the 
origen. 

Summarizing: as consequence of the renormalization 
procedure assumed in the Eq.(|58|), the local functions 
characterizing the local extensive subsystem are given by: 



Po (A,p, U) = p\np- ^pln f^r^ 



- p\nF (XU) 



e=^ (l-- zd,\nF(z)\ , i 
2A V 3 ^ >^z=(\v)^ 



So (e, /o; C/) = Ae - po (A, p, V) ■ 



(68) 



(69) 



(70) 



It is not difficult to predict some of the consequences 
of the above microscopic model. The energy threshold, 
U, is related with the Newtonian potential as follows: 



8 



U = —m(l> (r) , 



(71) 



and therefore, the vahie of U decreases from the inner 
regions of the system to the outer ones. In the inner re- 
gions the local subsystems exhibit the greater values of 
the parameter z, and therefore, their local distribution 
functions for the particles momenta are almost gaussian. 
However, at the outer regions, the local distribution func- 
tions diverge from the gaussian shape, becoming asymp- 
totically in a uniform distribution. It is seen that the 
gravity drives the behavior of the local distribution func- 
tion. 



IV. STRUCTURE EQUATIONS AND SYSTEM 
SCALING LAWS. 

From the first equilibrium relation of the Eqs.(|50|) it is 
deduced that the parameter A in the renormalized model 
is constant along all the system and it is equal to (3. 
The second equilibrium relation allows us to obtain the 
particle distribution, that is, the particle density: 



In /9 - 1 



2 Kh^pj 



I3m(j) + C, 



-InF 

that is: 

p = N{(3, fi) exp (-/3to(/) + C)F [(-/3m 
where the normalization constant N (/3, fi) is given by 

3 

'2m7r 



(72) 



7V(/3,/i) 



exp(-^- 1), 



(73) 



(It will be obviated the subindex s in the mean field de- 
scription). 

This last result allows us to specify the scaling param- 
eter X of the system scaling laws, Eqs.(p8|) and Eqs.(^0|), 
since from the Eq.([73|) the following relation is valid: 



and therefore: 



(74) 



(75) 



Taking into consideration the relations of the Eq.(28) 
the scaling exponent constants are given by: 



m = 2, rj = 



c=-i, x = 2, 



(76) 



and therefore, the energy scaling exponent r in the 
Eq.(|ll) is: 



(77) 



This result differs from the proposed by Vega & 
Sanchez in the refs. |3^j3^ . In these papers they stud- 
ied the present model but this time considering a box 
renormalization. They claimed correctly that the ther- 
modynamic limit for the astrophysical system must be 
invoked differently from the usual extensive systems. In 
the Eq.(l) of the ref. ||3^ they demand that the thermo- 
dynamic limit is performed when: 



N 



oo, E oo, L 



E N 
oo, — = const, — = const, 
N L 



(78) 



where L is the characteristic dimension of the box. How- 
ever, it is easy to see that this consideration leads to the 
non-proportionality of the Boltzmann's entropy with the 
number of particles (see Eq.(8) or the expression before 
the Eq.(12) in the ref. |3^), and therefore, ensemble in- 
equivalence. In spite of this error, all the results obtained 
by them are still correct, since they performed correctly 
the scaling invariant description with the introduction of 
the A'^independent parameters ^ and 77 defined as follows: 



EL 



(in the microcanonical ensemble), (79) 



Gm^N 
L 



(in the canonical ensemble). 



(80) 



which are also in agreement with the scaling laws consid- 
ered by us in the Eqs.(|^), (|l|), (^, and (|^). In their 
approach, they assumed incorrectly that the mechanical 
energy scales proportional to N. Our energy scaling law 
suggests that this model can not be applied for a num- 
ber arbitrarily large of particles: an upper bound for N 
appears for the validity of the non relativistic conditions 
(see in section VI). 

As it can be easily seen, in the Eq. (fz^) , when (p 0, 
the density vanishes, p 0. Thus, our model is renor- 
malized. In order to obtain the structure equations, the 
Eqs.(pl|), it is necessary to perform the following calcu- 
lation: 



d 



So (e,p; 



d 



pq (t>) 



z = (-f3m<p) 2 



Introducing the constant K: 
the dimensionless function <&: 



and the coordinate ^: 



-f3m(f). 



(81) 



(82) 



(83) 
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e = i^^r, (84) 
the structure equations are finally written as follows: 

A$ = -47rFi (C, $) , AC = -47ri^2 (C, $) , (85) 
where: 

Fi(C,$) =exp(C + $)i^($5^ , 



F2(C,$) = ^exp(C)$5 



(86) 



Due to the scaling transformation, the constant K de- 
pends on N as follows: 



K = N-^Kr, 



(87) 



where Ko is obtained substituting (3 by /3o in the Eq.(82) 
The coordinate ^ is A'^independent. 



V. NUMERICAL ANALYSIS 

In order to perform the numerical analysis of the above 
equations, it is convenient to introduce the following con- 
stant: 



2tt 



Po = (2^) 



3GW 



1 

2tt Gm? '' 



which are respectively the energy, density, and length 
characteristic units of the present analysis. The above 
consideration allows us to set Ti — G — m — \. Hearafter 
it will be imposed the condition iV = 1. 

In terms of the functions G and $, the particle number 
constrain is written as follows: 



3 ^00 

exp(-^-l):^ [ 47r^2^^Fi(C,$) = l, 

PS i 



(89) 



which establishes a functional relation between the ther- 
modynamic parameters Po and /i. This exigency can be 
rewritten analyzing the behavior of the functions G and $ 
in the asymptotic region ^ +00. The spherical solution 
of the Poisson's equation satisfy the following asymptotic 
relation: 



lim r^VrCt) (r) = GM, 



(90) 



where M is the total mass of the system. In the 
present problem, taking into consideration the charac- 
teristic units, the Eq.(p8|), the above relation is rewritten 
as follows: 



lim ^2^ci>(0 = -/3oi^| 

{—too Oc; 



(91) 



In the asymptotic region ^ — s- cx), the functions (f> and 
G vanish identically: 



hm <I>(0 = lim C(C) = 0, 

^^00 ^— *-CXD 



(92) 



due to the asymptotic behavior of the Green's function, 
Eq.(|o|), in their representations, Eq.(|l8|) and Eq.(|9|). 
The numerical solution of the Poisson's like equations 
system. Eg. (|85|) , is carried out imposing the boundary at 
the origin 1^ = 0. This is performed demanding at the 



origin: 



$(0) = $0 > 0, 
G(0) =Co. 



(93) 
(94) 



The value of Gq must be selected appropriately since 
it must vanish when ^ — > 00, and therefore, Gq depends 
on the parameter $o- This situation can be overcome 
redefining the problem as follows: firstly, displacing the 
function C (^): 



c(0 



(0. 



(95) 



where c ^ is the solution of the Poisson's like equation 
system with the following boundary conditions: 



$(0) =$0, 
c(0) =0, 



(96) 
(97) 



and the new coordinate ^ is related with ^ throughout 
the relation: 



£, = exp 



1 



c (00) 



Hereafter it is obviated the punctuation in ^. The 
Eq.(pl|) is rewritten as follows: 



h ($0) exp 



-c (co) 



PoKS = /3o" exp 



(a* + 1) 



where the function h (<&o) was introduced: 



.2 d 



/i($o) = - lim e'^'i>(0 



(99) 



(100) 



So far, there are three parameters involved in the so- 
lution of the problem, (3o, /i and $0. The above relation 
determines as a function of [3o and $o- An additional 
relation is needed to specify the functional dependency 
between Po and ^o- The solution to this situation is 
obtained throughout the equivalency between the sta- 
tistical ensembles when this equivalency held. From the 
Legendre's transformation is followed that the energy is 
obtained from the Planck's potential through the rela- 
tion: 



(101) 
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and therefore: 
dpo ($o) d 



(102) 



This is the second relation that we have been looking for. 

From the Eq.(|99|) it is deduced the functional depen- 
dency between the chemical potential n and the param- 
eters Po and 



H = -1 + - \n/3o - 21n/i ($0) - c (00) . 
Using the Eq. , the energy is expressed by: 



e(/3o,$o) 



1 

To 



3 _ /ii ($0) 
2 



(103) 



(104) 



where hi ($o) is defined as: 

+00 

hi ($0) = / 4^e'de$ 



F2(c,$)-K-i^i(c,$) 



(105) 



Similarly, the Planck's potential in the Eq.(|5J) is 
rewritten as follows: 



V (/3o, $0) = ln/3o + 2 In/j ($0) + ^J^^^: 

where h2 (^o) is defined by the expression: 
+00 



(106) 



($0) = y 4<2^cj^^ ^) 



(107) 



Introducing the fimction H ($0): 

i?(<i>o)=21nM*o) + ^)^: 
/i ($0) 

the Eq.( |l02|) is rewritten as follows: 



dln/3o($o) _ 5i7($o) 



3- 



fei (^0 

/i($o) 



(108) 



(109) 



(110) 



Of course, the validity of the Eq.( |l09| ) rests on the 
equivalency between the micro and canonical ensembles, 
which is satisfied wherever the parameter (3o is a decreas- 
ing function of the energy: 



Finally, the entropy of the system is given by: 
5 [e (/3„, $0)] = ^ + ^ In /3„ - if ($0) - 



dPo ($0) 
de ($0) 



< 0. 



(Ill) 



Introducing the function TZ {^q): 
<S>o 



7^($o) 



dHjs) 
ds 



Mf) 

h{s) 



ds, 



(112) 



the general solution of the Eq.( |l09[ ) is expressed as fol- 
lows: 



/?„($o) = Cexp[7^($o)], 



(113) 



where the integration constant C could be fixed as the 
unity without loss of generality. Using the function 
TZ ((f>o) and the function K. (<f>o): 



the condition of the Eq.(lll) is rewritten as 



a9o d^o 



(114) 



(115) 



The Eq.( |109| ) will be justified when the above relation 
is satisfied, otherwise, at first sight, the results obtained 
from this methodology has apparently no sense. Accord- 
ing to the results of the ref. , the above observation is 
correct: in the representation {£, N) the equivalency of 
the ensem ble d emand the validity of the condition given 
in the Eq.(lll). Nevertheless, it can he chosen another 
representation for the integrals of motion in which the 
equivalence between the ensemble takes place. For exam- 
ple, the canonical description could be performed intro- 
ducing the following representation for the integrals of 
motion: 

{£,N)^ {U^,N) , where = {£/N) , 

where Lp is an arbitrary piecewise monotonic function of 
e — £ /N at least two times differentiable with exception 
of the juncture boundary. As it can be easily seen, U,^ 
and £ possess the same scaling law in the thermodynamic 
limit N ^ 00, that is the reason why both representa- 
tions are adequate in the generalized canonical ensemble 
with probabilistic distribution function given by: 



1 



exp {-(3^ ■ 



(116) 



where {Pip, N) is the partition function in the present 
representation from which it is derived the corresponding 
Planck's potential: 



iPip,N) = -lnZ^(/3^,iV) 



(117) 



It is easy to see that in this case it is also valid the Leg- 
endre's transformation between the thermodynamic po- 
tentials when the thermodynamic limit is performed: 



(118) 
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with: 



where sb (e) is the Boltzmann's entropy per particle. In 
this representation change: 



(119) 



During the representation changes the canonical pa- 
rameter and the Planck' potential are transformed as 
follows: 



/3o ^ Pip — 



dip (e) 



de 



f3o 



V 



Vp=V 



(120) 

^W-^^Vo^, (121) 



\ de J \^ \ ^ 

in the way that the entropy remains unchanged: 

Si£,N)~^S{Up,N) = Si£,N). (122) 

In the new representation the equivalence of the nii- 
crocanonical and canonical ensembles takes place when 
it is satisfied the condition: 



^S{Up,N) = Pp—\nPp<0, (123) 

that is, the negative definition of the curvature tensor 
asociate to the new representation. This condition can 
be rewritten again as follows: 



de 



de 



( d^p {e 



V de 



f3o 



< 0. 



(124) 



Let us analyze a possible way for the selection of the 
function (p. The Eq.(124) can be rewritten as: 



de 



^ - j3 —\n( '^'^^^^ 
de ° de \ de 



de 



a{e), 



(125) 



where the function a (e) > 0, from the which is derived: 



dyje) 
de 



CPo (e) exp 



de 



(126) 



where C is a positive constant which could be set as unity. 
It is easy to see that a convenient choice for the function 
a (e), from the theoretical viewpoint, could be given by: 



-kie) 

a{e) = {PUe)=P^c 
k{e) 



if fc(e) < 0, 
const if k (e) = 0, 
if k (e) > 0, 



where it was introduced the nomenclature: 
dPo (e) 



kie) = 



de 



d' 

^SB (6) , 



(127) 



(128) 



de 



1 if k (e) < 0, 

= {PcexpiM if fc(e) = 0, 
/32 (e) if k (e) > 0. 



(129) 



As it can be seen, in this representation (p (e) is an 
piecewise monotonic function of e, and therefore, there is 
a bijective correspondence between e and (p{e). In this 
way it have been shown that it is always possible to chose 
a representation in which it is valid the correspondence 
among the ensembles allowing us to extend the appli- 
cability of the generalized canonical ensembles to those 
situations in which the ordinary heat capacity is negative. 

The reparametrization freedom of the microcanonical 
ensemble allows us to extend the validity of the mean field 
approach with an appropriate selection of the function 
ip. It can be easily shown that all the results of the 
canonical description in the different representations are 
only non-coincident when it is taken into account those 
observables involving second derivatives (or superior) of 
the thermodynamic potentials. That is the reason why 



it is considered that the Eq.(109) can be extended to 
those regions where there is an univocal dependence of 
the canonical parameter f3o with the energy e. 

In order to solve the Eq.(109) it is necessary to known 
the asymptotic behavior of the functions h ($0), hi ($0) 
and /i2 ('I'o) for small values of $0. This is done develop- 
ing the functions c and $ in a perturbative expansion of 
the parameter $0: 

$==¥+(5$ + ..., c = c + 5c+... (130) 

as well as the functions Fi (c, $) and F2 (c, <&): 

Fi (c,$) = Fi (c,¥) +(5Fi ((5c,(5$;c,$) + 

F2(c,$) = F2 (c, ¥)+ (5^2 ((5c, (5$; c, $) + ..., (131) 

where the problem given in the Eq.(|85|) is rewritten as 
follows: 

A$ -47rFi (c, $) , Ac = -AnFi (c, , (132) 
A ((5$) = -47r(5Fi {Sc, 5$; c, $) , (133) 



A ((5c) = ~4:TrSF2 {Sc, (5$; c, $) , . 



(134) 



Using the power expansion of the function F (z), 
Eq. (|64[) , it is very easy to obtain the following relations: 

Fi (c, $) = exp (c) :7^¥% F2 (c, $) = exp (c) 



6F1 ((5c, 5$; c, $) 



30F" 



exp (c) 



2—1 2—3 
5 3 



(135) 



(136) 
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5F2 ((5c, (5$; c, $) = exp ( 



_i 1 1 

<i>^(5c+-<f> -(5$ 



Similarly, the functions 5c and 5<I> are expressed as 
follows: 



In this asymptotic region, the structure equations (the 
Eq.(^5|)) become in the following problem: 

4 3 2 i 

A$ = -47rexp (c) — Ac = -47rexp (c) , 
Si/tt V^r 

(137) 



$ (0) = $0, $ (0) = and c (0) = 0, c' (0) = 0. (138) 

It is very easy to see that the above problem possesses a 
fractal characteristic and is quite similar to the polytropic 
model with polytropic index 7 = §, which allows us to 
express its general solution as follows: 



o^(4e), 5 = ^('i>Je), (139) 



$ = $ 



where the functions (z) and ip {z) are the solution of 
the problem: 

4 3 2 1 

A^v? = -47rexp (-0) -—^<f^, A^ip = -47rexp (V') —f='P^, 

(140) 

with boundary conditions: 

(0) = 1 and if' (0) = 0, V (0) = and V' (0) = 0. 

(141) 



It is interesting to point out that this polytropic in- 
dex characterizes an adiabatic process of the ideal gas of 
particles, which in our case is the evaporation of the sys- 
tem in the vacuum. However, this equation system is not 
equivalent to polytropic model due to the presence of the 
term of the gravity driving, the presence of the function c. 
It is easy to understand that the maximum effects of this 
term appear in the outer region of the system, in the halo. 
In this region the results of the polytropic model and the 
pseudo-polytropic model given by the Eq. ( |68|) are differ- 
ent. Moreover, this kind of behavior is also present for 
all the values of the parameter $0, and as a consequence 
of the renormahzation prescription assumed in Eq.(58), 
it leads to configurations of the systems characterized by 
an isothermal core with a quasi-polytropic halo. The ef- 
fect of the function c is showed in the FIG. 4. As it can be 
seen, the consideration of this term reduces the system 
size in comparison with the polytropic equation. 



FIG. 4. Comparison between the polytropic model with 
7=1 and the quasi-polytropic model presented in the present 
analysis. 



(5$ = $ 



where the functions ipi (z) and ipi (z) are obtained from 
the problem: 



i^zVi = -47rexp (ip) 



2 5 2 3 1 
5 3 



A,Vi 



-47rexp (V') 



(143) 



(144) 



with boundary conditions: 

ipi (0) = and (^'i (0) = 0, ipi (0) = and ip[ (0) = 0. 

(145) 

In this way it is obtained the asymptotic dependence 
of the functions h (s), hi (s) and /12 (s): 



h [s) = aisi + a2s3 + O {^s'^^ 
hi (s) = bisi + b2S^ + O (^s^^ 
h2 (s) = cisi + C2si + O {^si^ 
The calculation yields: 



ai = 0.744, 02 = -0.143, 

61 = 1.5ai = 1. 116, 62 = -0.242, 

ci = -0.622, C2 = 0.478. 



(146) 
(147) 
(148) 

(149) 



Thus, in the asymptotic region, the relation between 
the canonical parameter [3o and $0 is given by: 

/9o ($0) = $0 + + O ($^) , (150) 

where: 

g~ 0.025. (151) 
Finally, the energy in this limit is given by: 

e~0.055 + 0($o)- (152) 

VI. RESULTS AND DISCUSSIONS 

Let us begin the discussion analyzing the scaling laws 
of our system. In general way, the scaling laws determine 
the specific form of the thermodynamic formalism. The 
astrophysical systems exhibit exponential self-similarity 
scaling laws in the thermodynamic limit, which is the rea- 
son why the Boltzmann-Gibbs' Statistics is applicable to 
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the macroscopic description of this kind of systems in 
spite of that they are nonextensive. This analysis can 
be performed choosing an adequate selection of the rep- 
resentation of the integrals of motion. As it was shown 
in the section IV, the A^-body Newtonian sclf-gravitating 
gas exhibits self-similarity properties under the following 
thermodynamic limit: 

E 1 
N — !■ oo, E — > oo, L — > 0, — ^ — const., LN^ — const., 

(153) 

where L is a characteristic linear dimension of the system. 
This result differed from the obtained by Vega & Sanchez 
in refs. [||,|6). 

Many investigators do not pay so much attention to 
the analysis of the system scaling laws when it is per- 
formed its macroscopic description. Fortunately, most 
of the systems studied so far belong to the class of the 
pseudoextensive systems. When it is performed the A^- 
independent description, the results of the analysis do 
not depend on the scaling laws, and therefore, many of 
these results remain valid. However, the non or bad con- 
sideration of the scaling laws leads in many case to the 
trivial nonequivalence of the statistics ensembles. 

As example of the above affirmation it can be seen the 
anomalies presented in the dynamical study of the self- 
gravitating systems performed by Cerruti-Sola & Pcttini 
in the ref. [3^ ]. In that paper, they observed a weaken- 
ing of the system chaotic behavior with the increasing 
of the particles number N. It is very well-known the 
consequence of this fact on the ergodicity of the system: 
the chaotic dynamics provides the mixing property in 
the phase space necessary for obtaining the equilibrium. 
In that example the chaoticity time grows with N, and 
therefore, the systems could expend so much time to ar- 
rive to the equilibrium. Similar behavior has been seen in 
the dynamical study of the so called Hamiltonian Mean 
Field model (see for example in the ref. |Q). 

However, our analysis allows us to understand the ori- 
gin of this behavior, a least, for the self-gravitating gas. 
In that study it was considered that the energy is scaled 
proportional to A'^ during the realization thermodynamic 
limit, which is a false assumption. It is very easy to see 
that the dependency of the instability exponent Xh of 
the energy (Fig. 4. of the ref. p^) is corrected when it 
is considered the right scaling law of the energy: its pro- 
portionality to TVs. The anomalies disappear when it 
is assumed this scaling law for the energy. Anyway, in 
spite of they assumed a wrong scaling law for the energy, 
they obtained the correct dependency of the instability 
exponent A^f with the energy per particle: the power law 
£2. Our analysis suggests that the anomalies presented 
in the dynamical study of the Hamiltonian Mean Field 
model could possess the same origin. In future works we 
will consider this possibility . 

Some other consequences of this scaling laws are found 
analyzing the limit of applicability of the model. In the 



present work it was considered a classical gas of identi- 
cal particles with mass m. Taking into account the rest 
energy of the particles, the nonrelativistic limit is valid 
when the absolute value of the mechanical energy of the 
system is much smaller than its rest energy : 



eoe($o)Af5 ^mc^N, 



(154) 



where eg is the characteristic energy of the model, 
Eq. (psf) . However, this condition can not be satisfied for 
an arbitrary number of particles. In fact, when A'^ tends 
to No: 



A^n - 27r 



hc\ ^ 1 



(155) 



to which corresponds to a characteristic mass of Mq 

3 

1 



i2' 



(156) 



the model loses its validity. Everybody can recognize the 
fundamental constant of the stellar systems, which has 
much to do with the stability conditions of the stars (see 
in ref. Note that this constant appears as conse- 

quence of the energy scaling assumed, so that, it can not 
be obtained if a another scaling law had been adopted. 
A consequent analysis of these massive systems should 
be performed taking into account the relativisty effects. 

Finally, let us remember the well-known white dwarfs 
model based on the consideration of the Thomas-Fermi 
method to describe the state equation of the degenerate 
nonrelativistic electronic gas, whose pressure supports 
the hydrostatic equilibrium of the star. Using a sim- 
ple dimensional analysis, from this model can be easily 
derived the same scaling laws obtained by us analysing 
the necessary conditions for the equivalence of the statis- 
tical ensembles. This coincidence is not casual: for the 
nonrelativistic particles these scaling laws only depend 
on the dimension of the physical space. 

Let us discuss now the results of the numerical calcu- 
lations. In the FIG. 5 it is shown the general dependency 
of the functions h ($0), hi (<i>o) and /i2 ('i'o)- Observe the 
oscillatory character of these functions with the increas- 
ing of the parameter $0. This behavior is characteristic 
of the isothermal distribution. It can also be seen the 

quasi-polytropic asymptotic behavior for low values of 

3 

^0, the power law $q . 

FIG. 5. Behavior of the functions ^(<l>o), ^1 ("I'd) and 
h2 (*&o)- The oscillatory behavior is characteristic of the 
isothermal distribution. 



Throughout the Eq.( 109 ) it is obtained the dependence 
of the canonical parameter f3o, as well as the energy e 
from the parameter $o- Similarly, it is studied this de- 
pendence for the radio in which is contained the 80 % 
of the total mass of the system, i?80% j as well as for the 



14 



central density p^. All these dependencies are shown in 
the FIG. 6, 7, 8 and 9 respectively. 

FIG. 6. Canonical parameter /3 vs $0. Observe the persis- 
tence of the ascillatory behavior. 

FIG. 7. Scaling invariant energy, e vs $0. It is observe 
again the oscillatory behavior. 

FIG. 8. Radio at the 80% of the system total mass, -Rgo% 
vs $0. 

FIG. 9. Central density po vs $0. This function grows 
monotonically with the increasing of $0. This result allows 
us to understand that with the increasing of $0 the system 
develop a core-halo structure (a high dense core and a dilute 
halo) . 

Let US to comment briefly these results. At the first 
place it is noted the bounded character of the parameters 
/3o, e and i?80% for all the values of the parameter $0: all 
of them are contained in the following intervals: 

Po e (0, 2.97) ; e e (-0.26, 0.055) ; i?8o% e (1-02, 1.33) . 

(157) 

showing an oscillatory behavior for high values of $0 
around the values /3oo — 2.12, too — —0.18, and 
(i?80%)oo — l-ll- However, the central density grows 
monotonic with the increasing of $0i exhibiting an ex- 
ponential behavior for $0 > 9.45. It is very easy to un- 
derstand that our solution exhibits the same features of 
the Antonov problem solution . This is evident when 
analyzing the caloric curve: the dependency Po vs e, the 
FIG. 10: it is found again the very well-known spiral. Ad- 
ditionally, it is shown in the FIG. 11 the dependence of 
the radio in which is contained the 80% of the total mass 
of the system with the energy. 

FIG. 10. Caloric curve of the astrophysical system. It is 
obtained again the classical spiral. In general, the results of 
the present model exhibit the same features of the Antonov's 
problem. 

FIG. 11. System size vs Energy. The system becomes 
smaller with the energy decreasing. 

All these configurations contained between the points 
B and C can be described using the canonical description 
in the representation (5, N). Here the system exhibits a 
positive heat capacity. All the equilibrium configurations 
between the points C and D can not be accessed from 
the canonical description in the {£, N) representation. In 
this part of the spiral the systems exhibits a negative heat 
capacity. No equilibrium configurations exits for values of 



/?o greater than [3c = 2.97. Ordinarily this is referred as 
the isothermal catastrophe. On the other hand, no equi- 
librium configurations exits for energies outside the inter- 
val (—0.26,0.055). For energies greater than — 0.055 
the system is extreme diluted and can not be confined. 
For energies below of eo ~ —0.26, the systems collapse 
developing a core-halo structure. It is obtained again the 
gravothermal catastrophe. 

All those configurations contained between the points 
B to D correspond to equilibrium situations, the other 
points of the spiral correspond to unstable saddle points. 
This is evident when analyzing the thermodynamic po- 
tentials of the ensembles, the FIG. 12. and 13. For exam- 
ple, in the FIG. 12. it is represented the dependence of the 
Planck potential from the canonical parameter (3o. Here 
it is evident that in the canonical ensemble (using the 
(£, N) representation) all those configurations between 
the points B to C are stable: in such points it is min- 
imized this thermodynamic potential. Similarly, in the 
FIG. 13. it is shown that in the microcanonical ensem- 
ble all those configurations between the points B to D 
are stable: in such points it is maximized the entropy for 
each accessible value of the energy. 

FIG. 12. Planck potential vs canonical parameter /3. Here 
is evident that in the canonical ensemble are stable all those 
configurations between the points B to C. 

FIG. 13. In the microcanonical ensemble are stable all 
those configurations between the points B to D. 

It has to be pointed out that the Gibbs' argument, 
the equilibrium of a subsystem with a thermal bath, is 
not applicable to this situation, since it is based on the 
independence or weak correlation of the subsystem with 
the thermal bath. This is an invalid supposition for the 
nonextensive systems due to the long-range correlations 
existing in them as a consequence of the long-range in- 
teractions among the particles, which is the case that we 
are studying here. It is needed to remember that the use 
of generalized canonical description is supported by its 
equivalency with the microcanonical description in the 
thermodynamic limit, which is the description physically 
justified to this system. Moreover, no reasonable ther- 
mal bath exits for the astrophysical systems. Thus, the 
isothermal catastrophe is not a phenomenon with physi- 
cal relevance since it can never be obtained in nature: the 
consideration of a thermal bath in the astrophysical sys- 
tem is out of context. A different significance possesses 
the gravothermal catastrophe. The gravitational collapse 
is the main engine of structuration in astrophysics and 
it concerns almost all scales of the universe: the forma- 
tion of planetesimals in the solar nebula, the formation 
of stars, the fractal nature of the interstellar medium, 
the evolution of globular clusters and galaxies and the 
formation of galactic clusters in cosmology |]30| . 

Many authors claim that those regions characterized 
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by a negative heat capacity can not be accessed by the 
canonical description, and the only way to access to them 
is microcanonically . It should be remarked that this 
assumption is not strictly correct. The reason is easy to 
understand since the analysis of the equivalency between 
the canonical and microcanonical ensemble does not need 
to be performed using a global representation for the inte- 
grals of motions, the local representations are admissible 
too. The microcanonical ensemble is local reparametriza- 
tion invariant: the description of an arbitrary small re- 
gion of the integral of motion space could be performed 
using any representation of the integrals of motion: the 
physical observables in the microcanonical description do 
not depend on the representation (see in ref . ) . Thats 
is the reason why it is considered that the analysis of 
the equivalence of the microcanonical ensemble with the 
canonical one does not have to be limited to the use of 
only one representation for all the space (global analy- 
sis), but this analysis could be made using a different 
representation in an arbitrary small regions of the space 
(local analysis). In the previous section it was used this 



argument to justified the validity of the Eq.(109) for all 



the accessible equilibrium configurations of the system, 
showing in this way the importance of this geometrical 
characteristic of the probabilistic distribution functions 
in the macroscopic description of systems. 



VII. CONCLUDING REMARKS 

In the present paper it has been shown how the 
Boltzmann-Gibbs' Statistics can be improved in order to 
extend its applicability to the study of the macroscopic 
description of the A^-body self-gravitating gas, although 
this kind of systems are nonextensive. The Tsallis' Statis- 
tics is non applicable to this kind of system since this 
theory demands potential self-similarity scaling laws in 
the thermodynamic limit [ p^ , which is not the case that 
it is studying here: the astrophysical systems are pseu- 
doextensive ]l3[ , they exhibit an exponential scaling laws 
when N ^ oo. 

The most important result obtained in the present 
analysis is the specification of the thermodynamic limit 
for this kind of systems. Although the A^-independcnt 
picture of the traditional description of the astrophysi- 
cal systems remains inalterable, the bad consideration of 
the scaling laws leads to a trivial ensemble inequivalency, 
which was shown in the dynamical analysis performed in 
the ref. (37). 

The coincidences and connections of the results de- 
rived from this study with others results obtained in 
the past using both, thermodynamic and non thermo- 
dynamic methods, constitute important evidences of the 
validity of our considerations. 

Finally, it was shown the importance of the geometrical 
aspects of the probabilistic distribution functions in the 
macroscopic study of systems. Specifically, it was used 



these geometrical properties to extend the validity of the 
generalized canonical description to those situations in 
which the traditional description limited to the consid- 
eration of an special representation for the space of the 
integrals of motion of the system fails. These geometric 
aspects could be used to enhance the possibilities of the 
Montecarlo method based on the canonical exponential 
weight during the study of those hamiltonian systems 
presenting anomalies in their heat capacity. 



VIII. APPENDIX 

A. Derivation of the descomposition formula. 

The A^-body phase space Xn is an Cartesian external 
product of the generalized coordinates and momentum 
spaces of each particle, Q and V, which can be repre- 
sented in the following way: 



N 



(158) 



s=l 



Using the partition of the coordinates space Q in not 
aver lapped cells {c^}: 



Q =|J cfc, 



(159) 



the space Xn can be descomposed in not overlapped sub- 
Xn^\Jx^n\ (160) 



spaces X^ ' 



where 



N 



(161) 



Here, p^ (s) is a function which assigns each particle of 
the system to a determined cell. The index a denotes all 
the possible ways to perform thi s cor respondence. More- 
over, qs (fc) = Cfc. Using the Eq. (|l60| ), the A^-body phase 
space integration can be expressed as follows: 



1 

iV! 



dX 



N 



dX 



(162) 



A- : 



Due to the identity of particles, all those terms cor- 
responding to configurations with identical occupation 
number of particles at the cells are identical. Let {n^} 
be the occupation numbers of particles in the cells. For 
this this case, there are a total of: 



iV! 



(163) 
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identical terms in the sum of the Eq.(162). Let xjij^''^ 
be a subspace of with occupation numbers given by 
{ti/c}. Since Ajv is a Cartesian external product of spaces, 
its volume element can be factorized in the volume ele- 
ments of the different spaces. This can be convenientely 
done grouping all these particles beloging to the same 
cell.. In this case the A^-body phase space integration of 



[9; 

[10 

[11 

[12 



the bubopaco jYIJ^''^ can be doocompoood ao foUowo: 



M 



— ' /^iV — 



dx^ =n d 



X 



where: 



O 



i / dX„, if n 7^0, 
1, ifn = 0, 



(165) 



(k) 

and X„f, represents the Ti^-body phase space whose phys- 



ical Ispace of the particles is reduced to the cell Ck ■ Tak- 



[15 
[16 

[17: 

[18 
[19 



ing into account all the expossed above, it is easy to see 
that the iV-body phase space integration can be finally 
expressed as: 



Nl 



N 



Xn 



k Tit=0 



(166) 



where the 5'^'^^ (n) is related with the Kronekel delta func- 
tion as follows: 



<5(^) in) = Son. 



(167) 



The presence of this function in the right hand of the 
Eq.( |166[ ), assures that the number of particles remain 
fixed and equals N. 
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